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Abstract 

For m number of bosons, carrying spin (s = |) degree of freedom, in VL number of single 
particle orbitals, each doubly degenerate, we introduce and analyze embedded Gaussian orthogonal 
ensemble of random matrices generated by random two-body interactions that are spin (S) scalar 
[BEGOE(2)-s]. Embedding algebra for the BEGOE(2)-s ensemble and also for BEGOE(l+2)-s 
that includes the mean-field one-body part is C7 (20) D U (Q) <g> SU(2) with SU{2) generating spin. 
A method for constructing the ensembles in fixed- (m, 5) spaces has been developed. Numerical 
calculations show that for BEG0E(2)-s, the fixed-(m, S) density of states is close to Gaussian and 
level fluctuations follow GOE in the dense limit. For BEG0E(l+2)-s, generically there is Poisson 
to GOE transition in level fluctuations as the interaction strength (measured in the units of the 
average spacing of the single particle levels defining the mean-field) is increased. The interaction 
strength needed for the onset of the transition is found to decrease with increasing S. Propagation 
formulas for the fixed- (m, S) space energy centroids and spectral variances are derived for a general 
one plus two-body Hamiltonian preserving spin. Derived also is the formula for the variance 
propagator for the fixed- (m, S) ensemble averaged spectral variances. Using these, covariances 
in energy centroids and spectral variances are analyzed. Variance propagator clearly shows, by 
applying the Jacquod and Stone prescription, that the BEG0E(2)-s ensemble generates ground 
states with spin S = S max . This is further corroborated by analyzing the structure of the ground 
states in the presence of the exchange interaction S 2 in BEG0E(l+2)-s. Natural spin ordering 
(Smax, S max — 1, S max — 2, . . ., or i) is also observed with random interactions. Going beyond 
these, we also introduce pairing symmetry in the space defined by BEG0E(2)-s. Expectation 
values of the pairing Hamiltonian show that random interactions exhibit pairing correlations in the 
ground state region. 

PACS numbers: 05.30.Jp, 05.45.Mt, 03.65.Aa, 03.75.Mn 
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I. INTRODUCTION 



Random matrix theory has been established to be one of the central themes for chaotic 
quantum systems [1]. The classical Gaussian orthogonal (GOE), unitary (GUE) and sym- 
plectic (GSE) ensembles of random matrices are ensembles of multi-body, not two-body 
interactions. However, finite interacting quantum systems, such as nuclei, atoms, quan- 
tum dots, small metallic grains and interacting spin systems modeling quantum computing 
core and BEC, are governed largely by two-body interactions and hence, it is important to 
consider ensembles generated by random two-body interactions. These ensembles are de- 
fined by representing the two-particle Hamiltonian by one of the classical ensembles (GOE, 
GUE, GSE) and then the many particle (to > 2) Hamiltonian is generated by exploiting 
the direct product structure of the m-particle Hilbert spaces. As a random matrix ensemble 
in the two-particle spaces is embedded in the many particle Hamiltonian, these ensembles 
are generically called embedded ensembles (EE) and with GOE embedding, they will be 
EGOE's. A wide variety of EE for fermions have been introduced in literature [2-7]. 

Simplest of EE is the embedded Gaussian orthogonal ensemble for spinless fermion sys- 
tems generated by random two-body interactions, denoted by EGOE(2). In general, it is also 
possible to define EGOE(/c) ensembles generated by /c-body (k < m) interactions [4]. It is 
useful to mention that many diversified methods like numerical Monte-Carlo methods, binary 
correlation approximation, trace propagation, group theory, supersymmetry and perturba- 
tion theory are used to derive generic properties of EE [4, 8-11]. Some of the generic results 
for EGOE(/c) are as follows: (i) eigenvalue density exhibits, with increasing to, transition 
from semicircle to Gaussian with to = 2k being the transition point [8]; (ii) numerical studies 
have shown that the level and strength fluctuations follow GOE [12]; (iii) there is average- 
fluctuation separation with increasing to [12, 13]; (iv) ensemble averaged transition strength 
densities follow bivariate Gaussian form and consequently, transition strength sums will be 
close to a ratio of two Gaussians [14]; and (v) there will be non-zero correlations between 
states with different particle numbers [15, 16]. Besides two-body interactions, Hamiltonians 
for realistic systems contain a mean-field part [defined by non-degenerate single particle (sp) 
levels] and then the ensemble is denoted by EE(l+2). This ensemble exhibits, with increas- 
ing strength A of the interaction (A is in units of the average sp level spacing), three chaos 
markers defining transitions in level fluctuations (A c ), strength functions (Xp) and entropy 
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(Xd) respectively. Generic results of EE(2) are valid for EE(l+2) in the strong coupling 
regime (A >> A^). Realistic systems also preserve various symmetries. For example, spin S 
is a good quantum number for atoms and quantum dots, angular momentum J and parity 
7r are good quantum numbers for nuclei and so on. Therefore, it is more appropriate to 
study EE with good symmetries. A simple but non-trivial extension of EE is to consider 
finite quantum systems with spin (s = |) and then the ensembles are called EE(l+2)-s. 
For finite interacting fermion systems, EG0E(l+2)-s has been studied in detail using a 
mixture of numerical and analytical techniques [6, 17, 18] and EGUE(2)-s using Wigner- 
Racah algebra [10]. Moreover, EGUE(2) with spin-isospin SU(A) symmetry, comparatively 
complicated than the random matrix models analyzed before, has been recently introduced 
and analyzed in some detail [11]. More importantly, there are now several applications of 
EE to mesoscopic systems [19-22], quantum information science [23] and in investigating 
thermalization in finite quantum systems [24-27]. Unlike for fermion systems, there are only 
a few EE investigations for finite interacting boson systems [28-32]; the corresponding EE 
are called BEE (B stands for bosons). Briefly, these studies are as follows. 

Firstly, it is important to mention that, unlike fermion systems, for interacting spinless 
boson systems with m bosons in N sp orbitals, dense limit defined by m — » oo, iV — > oo and 
m/N — > oo is also possible as m can be greater than N for bosons. It is now well understood 
that BEGOE(2) [also BEGUE(2)] generates in the dense limit, eigenvalue density close to 
a Gaussian [28, 33]. Also the ergodic property is found to be valid in the dense limit with 
sufficiently large N [31]; there are deviations for small N [29]. Similarly, for BEGOE(l+2), 
as the strength A of the two-body interaction increases, there is Poisson to GOE transition 
in level fluctuations at A = A c [31] and with further increase in A, there is Breit-Wigner 
to Gaussian transition in strength functions [32]. For BEGUE(A;), exact analytical results 
for the lowest two moments of the two-point function have been derived by Agasa et al 
[30]. Level fluctuations and wavefunction structure in interacting boson systems are also 
studied using interacting boson models of atomic nuclei [34-36] and a symmetrized two 
coupled rotors model [37, 38]. In addition, using random interactions in interacting boson 
models, there are several studies on the generation of regular structures in boson systems 
with random interactions [39-42]. Finally, there are also studies on thermalization in finite 
quantum systems using boson systems [24-27]. 

Going beyond the embedded ensembles for spinless boson systems, our purpose in this 
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paper is to introduce and analyze spectral properties of embedded Gaussian orthogonal 
ensemble of random matrices for boson systems with spin degree of freedom [BEG0E(2)-s 
and also BEG0E(l+2)-s] and for Hamiltonians that conserve the total spin of the m-boson 
systems. Here the spin is, for example, as the F-spin in the proton-neutron interacting 
boson model (pnlBM) of atomic nuclei [43]. Just as the earlier embedded ensemble studies 
[29-32], a major motivation for the study undertaken in the present paper is its possible 
applications to ultracold atoms. There are several studies of the properties of a mixture of 
two species of atoms which correspond to pseudospin-| bosons (i.e. two-component boson 
systems) with m s = ±| distinguishing the two species; see for example [44, 45]. However, 
the Hamiltonians appropriate for these studies do not conserve the total spin (as the system 
does not have true |-spins) and therefore, the model study presented in the present paper 
will not be directly applicable to these systems in understanding their statistical properties. 
Nevertheless, the BEGOE(l+2)-s with spin-| bosons is a simple yet non-trivial extension 
of the spinless BEE. This ensemble is useful in obtaining several physical conclusions, like 
spin dependence of the order to chaos transition marker in level fluctuations, the spin of 
the ground state (gs), the spin ordering of excited states and pairing correlations in the 
gs region generated by random interactions, that explicitly require inclusion of spin degree 
of freedom (these are discussed in Sections III, V and VI). It should be emphasized that 
the present paper opens a new direction in defining and analyzing embedded ensembles for 
boson systems with symmetries. It is important to mention here that there are now many 
studies of spinor BEC using Hamiltonians conserving the total spin with the bosons carrying 
s = 1 (also higher) degree of freedom [46, 47]. Extensions of BEGOE(l+2)-s with s = \ 
to boson ensembles with integer spin s = 1 (or higher) is for future. Now we will give a 
preview. 

In Section II, introduced is the new embedded ensemble BEGOE(2)-s [and also BEGOE(l+2)- 
s] for a system of m bosons in Q number of sp orbitals that are doubly degenerate with total 
spin S being a good symmetry. A method for the numerical construction of this ensemble 
in fixed- (m,S) spaces is described. Numerical results for the ensemble averaged eigenvalue 
density, nearest neighbor spacing distribution and the long-range rigidity measure A 3 are 
presented in Section III. Propagation formulas for fixed-(m, S) energy centroids and spectral 
variances for general one plus two-body Hamiltonians that preserve S are given in Section 
IV. Here, given also is the analytical formula for the ensemble averaged fixed-(m, S) spectral 
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variances. Using these, studied are covariances in energy centroids and spectral variances 
generated by BEGOE(2)-s ensemble between states with different (to, £). Section V gives 
results for the preponderance of maximum £-spin ground states and natural spin order gen- 
erated by random interactions. Here, exchange interaction is added to the BEGOE(l+2)-s 
Hamiltonian. Pairing in BEGOE(2)-s is introduced in Section VI and presented also are 
some numerical results for pairing correlations. Finally, Section VII gives conclusions and 
future outlook. 

II. DEFINITION AND CONSTRUCTION OF BEGOE(l+2)-s 

Let us consider a system of to (to > 2) bosons distributed in Vt number of sp orbitals 
each with spin s = \. Then the number of sp states is N = 2Vt. The sp states are denoted 
by \i,m s = ±| ) with i = 1,2, ... ,Q and the two particle symmetric states are denoted by 
\(ij)s,m s ) with s = or 1. It is important to note that for EGOE(l+2)-s, the embedding 
algebra is U(2Q) D U(Q) <g> £77(2) with £77(2) generating spin; see Sections V and VI ahead. 
The dimensionalities of the two-particle spaces with s = and s = 1 are Q(Q — l)/2 and 
Q(Q + l)/2 respectively. For one plus two-body Hamiltonians preserving to particle spin £, 
the one-body Hamiltonian is h(l) = Yl?=i e i n i where the orbitals i are doubly degenerate, 
are number operators and are sp energies (it is in principle possible to consider h(l) with 
off-diagonal energies e^). The two-body Hamiltonian V{2) preserving to particle spin £ is 
defined by the symmetrized two-body matrix elements Vfj kl = (^(kl)s,m s \ V(2) \ {ij)s,m^j 
with s = 0, 1 and they are independent of the m s quantum number; note that for s = 0, 
only i 7^ j and k ^ I matrix elements exist. Thus V(2) = U s=0 (2) + V s=1 (2) and the sum 
here is a direct sum. The BEGOE(2)-s ensemble for a given (to, £) system is generated 
by first defining the two parts of the two-body Hamiltonian to be independent GOEs in 
the two-particle spaces [one for V s=0 (2) and other for U s=1 (2)], with the matrix elements 
variances being unity (except for diagonal matrix elements whose variance is two). Now the 
U(2) ensemble defined by {U(2)} = {V s=0 {2)} + {V s=1 (2)} is propagated to the (to, £)- 
spaces by using the geometry (direct product structure) of the m-particle spaces; here { } 
denotes ensemble. By adding the h(l) part, the BEGOE(l+2)-s is defined by the operator 

WBEGOE(l+2)-s = M 1 ) + A o {U s= °(2)} + X 1 {V s =\2)} . (1) 
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Here A and Ai are the strengths of the s = and s = 1 parts of V(2) respectively. The 
mean-field one-body Hamiltonian h(l) in Eq. (1) is defined by sp energies €j with average 
spacing A. Without loss of generality, we put A = 1 so that A and Ai are in the units of 
A. In the present paper, we choose sp energies q = i + 1/i as in our previous studies of 
fermion systems [6]. In principle, many other choices for the sp energies are possible. Thus 
BEGOE(l+2)-s is defined by the five parameters to, S, A , Ai). The H matrix dimension 
d(Q, m, S) for a given (to, S) is 

Mom** (2S+l) fn + m/2 + S-l\ (Q + m/2-S-2\ 

d ^ m > S) = ln^T){ m/2 + S + l ){ m/2-S J' (2) 

and they satisfy the sum rule J2s ( 2S + X ) d(n,m,S) = { N+ ^)- For example: (i) 
d(4, 10, S) = 196, 540, 750, 770, 594 and 286 for spins 5 = 0-5; (ii) d(4, 11, S) = 504, 900, 
1100, 1056, 780 and 364 for S = 1/2- 11/2; (iii) d(5, 10, S) = 1176, 3150, 4125, 3850, 2574 
and 1001 for S = - 5; (iv) d(6,12,S) = 13860, 37422, 50050, 49049, 36855, 20020 and 
6188 for S = - 6; and (v) d(6, 16, S) = 70785, 198198, 286650, 321048, 299880, 235620, 
151164, 72675 and 20349 for S = - 8. 

Given and V^ M , the many particle Hamiltonian matrix for a given (to, 5) can be 
constructed using the representation (Ms is the quantum number) and for spin 
projection the S 2 operator is used as it was done for fermion systems in [17]. Alternatively, 
it is possible to construct the H matrix directly in a good S basis using angular-momentum 
algebra as it was done for fermion systems in [48]. We have employed the Ms representation 
for constructing the H matrices with M s = M™ ra = for even to and M s = M™ n = \ 
for odd to and they will contain states with all S values. The dimension of this basis space 
is P(fi,TO,M™ n ) = J2s d(tt,m,S). For example, Z>(4, 10, 0) = 3136, £>(4,11,|) = 4704, 
£>(5, 10, 0) = 15876, P(6, 12, 0) = 213444 and P(6, 16, 0) = 1656369. 

To construct the many particle Hamiltonian matrix for a given (to, S), first the sp states 
\i,m a = are arranged in such a way that the first states have m s = | and the 
remaining f2 states have m s = —\ so that the sp states are |r) = \i — r, m s — |) for 
r < VL and |r } = \i = r — il, m s = — |) for r > il. Using the direct product structure of the 
many-particle states, the m-particle configurations m, in occupation number representation, 
are 

N=2fl \ 

m = ] [ m r ) = \m 1 ,m 2 , . . . ,m n ,m n+1 ,m n+2 , ■ ■ ■ ,m 2 n) , (3) 



where m r > with Ylr=i m r — m an d Ms = \ Y^l=\ m r ~ Ylt^Ln+i m fJ • To proceed fur- 
ther, the (l+2)-body Hamiltonian defined by q and V^ 0, is converted into the \i,m s = ±| ) 
basis. Then the sp energies d i with % — 1, 2, . . . , N are d { = e' i+n = ti for i < Q. Similarly, 
Vf jkl are changed to V imujmjtkmkMl = (im^jm^ V(2)\ km k ,lmi) using, 



^2'-? 2'^2'^2 ^ 2'^ 2'^ 2'^ 2^ ^ijkl ? 



(4) 



with all the other matrix elements being zero except for the symmetries, 

Using (e' r , V imu j mj ^m k ,imi)'' s i construction of the m-particle H matrix in the basis defined by 
Eq. (3) reduces to the problem of BEGOE(l+2) for spinless boson systems and hence Eq. 
(4) of [28] will give the formulas for the non-zero matrix elements. For completeness, these 
formulas are given in Appendix A. Now diagonalizing the S 2 matrix in the basis defined by 
Eq. (3) will give the unitary transformation required to change the H matrix in Ms basis 
into good S basis. Following this method, we have numerically constructed BEGOE(l+2)- 
s in many examples and analyzed various spectral properties generated by this ensemble. 
In addition, we have also derived some analytical results as discussed ahead in Sections 
IV and VI. These results are also used to validate the BEGOE(l+2)-s numerical code we 
have developed. In addition, we have also verified the code by directly programming the 
operations that give Eq. (A3). In this paper, we deal mainly with BEGOE(2)-s and the 
focus is on the dense limit defined by m - > oo, f2 — > oo, m/Q — > oo and S is fixed. Now we 
will discuss these results. 



III. NUMERICAL RESULTS FOR EIGENVALUE DENSITY AND LEVEL FLUC- 
TUATIONS IN THE DENSE LIMIT 

We begin with the ensemble averaged fixed-(m, S) eigenvalue density p m ' s (E), the one- 
point function for eigenvalues. First we present the results for BEGOE(2)-s ensemble defined 
by h{l) = in Eq. (1) and then the Hamiltonian operator is, 

WBEGOE(2)-s = A o {V s=0 (2)} + Ai {V s=1 (2)} . (6) 
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We have considered a 500 member BEGOE(2)-s ensemble with f2 = 4 and m — 10 and 
similarly a 100 member ensemble with Q = 4 and m = 11. Here and in all other numerical 
results presented in the paper, we use A = Ai = A. In the construction of the ensemble av- 
eraged eigenvalue densities, the spectra of each member of the ensemble is first zero centered 
and scaled to unit width (therefore the densities are independent of the A parameter). The 
eigenvalues are then denoted by E. Given the fixed- (to, S) energy centroids E c (m,S) and 
spectral widths a(m, S), E = [E — E c (m, S)\/o(m, S). Then the histograms for the density 
are generated by combining the eigenvalues E from all the members of the ensemble. Re- 
sults are shown in Fig. 1 for a few selected S values. The calculations have been carried 
out for all S values (the results for other S values are close to those given in the figure) and 
also for many other BEGOE(2)-s examples. It is clearly seen that the eigenvalue densities 
are close to Gaussian (denoted by Q below) with the ensemble averaged skewness (71) and 
excess (72) being very small; I71 1 ~ 0, I72 1 ~ 0.1 — 0.27. The agreements with Edgeworth 
(ED) corrected Gaussians are excellent. The ED form that includes 71 and 72 corrections is 
given by 



Ved(E) = rjg(E) { 1 + 



Vg(E) = -4=exp 
v 2ir 




(7) 



Here, He are Hermite polynomials: He^x) = x 3 — 3x, He^(x) = x A — 6x 2 + 3 and He%(x) = 
x 6 - 15x 4 + 45a; 2 - 15. 

For the analysis of level fluctuations (equivalent to studying the two-point function for 
the eigenvalues), each spectrum in the ensemble is unfolded using a sixth order polyno- 
mial correction to the Gaussian and then the smoothed density is r](E) = r]g(E){l + 
Ylf>3(0)' 1 S^He^(E)} with Co = 6 [13, 28]. The parameters are determined by min- 
imizing A 2 = J2 d iS ,m,S) [ F (Ei) - F(E)} 2 . The distribution function F(E) = rj(x)dx and 



similarly F(E) is defined. We require that the continuous function F(E) passes through 
the mid-points of the jumps in the discrete F(E) and therefore, F(Ei) = (i — 1/2). The 
ensemble averaged Arms is ~ 3 for ( = 3, ~ 1 for Co = 4 and ~ 0.8 for Co = 6 with some 
variation with respect to S. As Arms ~ 0.88 for GOE, this implies GOE fluctuations set in 
when we add 6th order corrections to the asymptotic Gaussian density. Using the unfolded 
energy levels of all the members of the BEGOE(2)-s ensemble, the nearest neighbor spac- 
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ing distribution (NNSD) that gives information about level repulsion and the Dyson-Mehta 
A 3 (L) statistic that gives information about spectral rigidity are studied. Results for the 
same systems used in Fig. 1 are shown in Fig. 2 for S = 2 and 5 (for other spins, the results 
are similar). In the calculations, middle 80% of the eigenvalues from each member are em- 
ployed. It is clearly seen from the figures that the NNSD are close to GOE (Wigner) form 
and the widths of the NNSD are ~ 0.288 (GOE value is ~ 0.272). The Aj^(L) values show 
some departures from GOE for L > 30 for S = S max and this could be because the matrix 
dimensions are small for S = S max in our examples (also the systems considered are not 
strictly in the dense limit and numerical examples with much larger m and Q with m » Q 
are currently not feasible). It is useful to add that S = S max states are important for boson 
systems with random interactions as discussed in Sections IV- VI ahead. In conclusion, sixth 
order unfolding removes essentially all the secular behavior and then the fluctuations follow 
closely GOE. This is similar to the result known before for spinless boson systems [13, 31]. 

Going beyond BEGOE(2)-s, calculations are also carried out for BEGOE(l+2)-s systems 
using Eq. (1) with Ao = Ai = A. We have verified the Gaussian behavior for the eigenvalue 
density for BEGOE(l+2)-s; an example is shown in Fig. 3a. This result is essentially 
independent of A. In addition, we have also verified that BEGOE(l+2)-s also generates 
level fluctuations close to GOE for A 0.1 for f2 = 4 and m — 10, 11 systems. Figure 3 
shows an example with A = 0.1. Going beyond this, in Fig. 4, we show the NNSD results, 
for a 100 member BEGOE(l+2)-s ensemble with Q = 4, m — 10 and total spins S — 0, 2 
and 5, for A varying from 0.01 to 0.1 to demonstrate that as A increases from zero, there is 
generically Poisson to GOE transition. A similar study has been reported in [6] for fermion 
systems. As discussed there, for very small A, the NNSD will be Poisson (as we use sp 
energies to be e; = % + 1/i, the A = limit will not give strictly a Poisson). Moreover, as 
discussed in detail in [6], the variance of the NNSD can be written in terms of a parameter A 
(A is a parameter in a 2x2 random matrix model that generates Poisson to GOE transition) 
with A = giving Poisson, A ^ 1 GOE and A = 0.3 the transition point A c that marks 
the onset of GOE fluctuations. We show in Fig. 4, for each A, the deduced value of A 
from the variance of the NNSD (Fig. 2 gives results for A — > oo for S = 2 and 5). As 
seen from the Fig. 4, A c = 0.039, 0.0315, 0.0275 for S = 0, 2 and 5 respectively. Thus A c 
decreases with increasing spin S and this is opposite to the situation for fermion systems. 
For a fixed value, as discussed in [6], the A c is inversely proportional to K, where K is the 
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number of many-particle states [defined by h(l)] that are directly coupled by the two-body 
interaction. For fermion systems, K is proportional to the variance propagator but not for 
boson systems as discussed in [31]. At present, for BEGOE(l+2)-s we don't have a formula 
for K. However, if we use the variance propagator Q(Q, m, S) for the boson systems [see Eq. 
(14) and Fig. 5 ahead], then qualitatively we understand the decrease in A c with increasing 
spin. 

Finally, it is well known that the Gaussian form for the eigenvalue density is generic for 
embedded ensembles of spinless fermion [5] and boson [31, 33] systems. In addition, ensemble 
averaged fixed-(m, S) eigenvalue densities for the fermion EGOE(l+2)-s are shown to take 
Gaussian form in [6, 17]. Hence, from the results shown in Figs. 1 and 3a, it is plausible 
to conclude that the Gaussian form is generic for EE (both bosonic and fermionic) with 
good quantum numbers. With the eigenvalue density being close to Gaussian, it is useful to 
derive formulas for the energy centroids and ensemble averaged spectral variances. These in 
turn, as discussed ahead, will also allow us to study the lowest two moments of the two-point 
function. From now on, we will drop the 'hat'over the operators H, h(l) and V(2). 



IV. ENERGY CENTROIDS, SPECTRAL VARIANCES AND ENSEMBLE AVER- 
AGED SPECTRAL VARIANCES AND COVARIANCES 



A. Propagation formulas for energy centroids and spectral variances 



Given a general (l+2)-body Hamiltonian H = h(l) + V(2), which is a typical member of 
BEGOE(l+2)-s, the energy centroids will be polynomials in the number operator and the 
S 2 operator. As H is of maximum body rank 2, the polynomial form for the energy centroids 
is (H) m,s = E c (m, S) = a + aim + a 2 m 2 + a 3 S(S + 1). Solving for the a's in terms of the 
centroids in one and two particle spaces, the propagation formula for the energy centroids 
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is, 



(H) m ' s = E c (m, S) = \(h(l))^] m + X ((^=°(2))) 2 ' J^^y 



+ A! ((V- =1 (2)» 



4Q{Q + 1) ' 
P°(m,S) = [m(m + 2)-AS(S + l)] , 

P\m,S) = [3m(m-2)+AS(S+l)} , ^ 

((v^(2))) 2 >° = E 1 «^ 1 (2)» 2 ' 1 = E V*g . 

i<j i<j 

For the energy centroid of a two-body Hamiltonian [member of a BEG0E(2)-s], the h(l) 
part in Eq. (8) will be absent. 

Just as for the energy centroids, polynomial form for the spectral variances 

cr H=h(i)+v(2)( m , s ) = (H 2 ) m,S - [E c (m, S)] 2 

is J2p=o a p mP + Yl 2 q =o b q m q S(S + 1) + co[S(S + l)] 2 . It is well known that the propagation 
formulas for fermion systems will give the formulas for the corresponding boson systems 
by applying Q — > —Q transformation [9, 33, 49-51]. Applying this transformation to the 
propagation equation for the spectral variances for fermion systems with spin given by Eq. 
(8) of [17], we obtain the propagation equation for &H=h(i)+v(2)( rn 'i * n terms of inputs that 
contain the single particle energies defining h(l) and the two particle matrix elements V^ kl . 
The final result is, 

4=MD+v( 2 )K S) = (HT' S ~ [E c (m, S)] 2 
(n-2)mm* + 20 (S 2 ) 



(Q-1)Q(Q + 1) ^ 

m*P°(m,S) *, n x 
7 -1 YOYO. + 1^1 2^ e ^A u ) 



„,*po 

+ 2(0~- 1)0(0 + 1) 

(O- 2)m*P\m, S) + 80(m - 1) (S 2 ) 
+ 2(0-l)0(0 + l)(0 + 2) ^ 6i m( } 

+P^=°(m,S) ^.(0) + P-^KS) E^W 

+ 4(o - 1)0(0 + i)(o + 2) 2rj Al ^ )Al ^ l} 
+p»= 2 > s =%m,s) ((y»= 2 > s=0 f) 2, \p u = 2 > s = l (m,s) ((V u = 2 > s=l ) 



(9) 
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The propagators P u ' s, s, which are used later, are 
P^=%m, S) = [(m + 2)m72 " ( ^ )] P ° (m ' S) 



(10) 



8(n-2)(n- i)n(n + 1) 

s= 8^(m - + 2m - 4) (g) + - 2)P 2 (m, ^(m, g) 

lm ' j 8(fi-l)fi(fi + l)(fi + 2) 2 

P^= 2 ' s=0 (m, 5) = [m*(m* - 1) - (S 2 )] P°(m, S)/[8Q{Q + 1)] , 

P^= 2 > s=1 (m, S) = {[(S 2 )] 2 (3fi 2 + 7tt + 6)/2 + 3m(m - 2)m*(m* + 1) x 

(ft - - 2)/8 + [(5 2 ) /2] [(5ft + 3)(ft - 2)mm* + ft(ft - l)(ft + l)(ft - 6)] }/ 

[(ft- i)ft(ft + 2)(ft + 3)] ; 

P 2 (m, 5) = 3(m - 2)m*/2 + (S 2 ) , m* = ft + m/2 , (S 2 ) = S(S + 1) . 
The inputs in Eq. (9) are given by, 



K,i( s ) = E % (! + " E (! + *«) > 

i k,l 

M«) = E V / ( 1 + ^)( 1 + ^)^ for i*j, 

k 

v m' s = v m ~ [(^( 2 )) 2 ' s + (A M ( 8 ) + A J - J ( 8 ))(n-2 (-l)')- 1 ] , 

= - (n - st-irr 1 v ' (l • '^Hi • ^)K, for < ^ j , 

v m' S = V ijki fOT all other cases . 

Eqs. (8) and (9) can be applied to individual members of the BEGOE(l+2) ensemble. On 
the other hand, it is possible to use these to obtain ensemble averaged spectral variances 
and ensemble averaged covariances in energy centroids just as it was done before for fermion 
systems [6]. Now we will consider these. 



(11) 



B. Ensemble averaged spectral variances for BEGOE(2)-s 

In this subsection, we restrict to H = V(2), i.e. BEGOE(2)-s and consider BEGOE(l+2)- 
s at the end. 

For the ensemble averaged spectral variances generated by H, only the fourth, fifth, 
seventh and eighth terms in Eq. (9) will contribute. Evaluating the ensemble averages of 
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the inputs in these four terms, we obtain, 



E Aj,-(o) = Ag(n-i)(n-2)(n + 2), 



E A?,.(l) = A?(n-l)(n + 2) 2 , 

id (12) 

Note that these inputs follow from the results for EGOE(2)-s for fermions given in [6] by 
interchanging s = with s — 1. Now the final expression for the ensemble averaged variances 
is 

a'tim.S) = ^Aj(!i-l)(!i-(-l)«2)(!i + 2)P"- 1 '«(m,S) 

+ ?"-y + l) ii+2) p-^v, s) (i3> 



In most of the numerical calculations, we employ Ao = Ai = A and then ajj(m, S) takes the 
form, 

a 2 H (m, S) Ao ^> =A A 2 Q(fi, m, S) . (14) 

Expression for the variance propagator Q(Q,,m,S) follows easily from Eqs. (8), (10) and 
(13). In Fig. 5, we show a plot of Q(fl, m, S)/Q(Q, m, S rnax ) vs S/S max for various Q and 
m values. It is clearly seen that the propagator value increases as spin increases and this 
is just opposite to the result for fermion systems [6]. An important consequence of this is 
BEGOE(2)-s gives ground states with S = S max [for fermion EGOE(2)-s, the ground states 
with random interactions have 5 = 0]. This result follows from the Jacquod and Stone [21] 
criterion and according to this (from the assumption of Gaussian form for the eigenvalue 



densities), the gs energy E gs is given by E gs oc — y <J 2 H {m, S). 

Before proceeding further, let us remark that for the BEGOE(l+2)-s Hamiltonian {if} 



h(l) + {V(2)}, assuming that h(l) is fixed, we have a 2 H = + <J v(2) • ^ ne ^ rs ^ term cr^^ is 
given by the first term of Eq. (9) and the second term is given by Eq. (13). In the situation 



h(l) is represented by an ensemble independent of {V(2)}, we have to replace by a 2 .^ 



in a 2 H . 
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C. Ensemble averaged covariances in energy centroids and spectral variances for 
BEGOE(2)-s 



Normalized covariances in energy centroids and spectral variances are defined by 

(H) m ' s (H) m '' s ' - \(H) m > s ) \(H) m '> s '\ 
Zu(m,S :m',S') = 1 J 1 J 



£ 22 (m, S : to', S') 



\h T ' s } {(#T'' 5 '} 

(Hi) m ' S {H 2 ) m '' S ' - j(# 2 > m ' 5 } {(# 2 > m '' 5 '} 



(15) 



These define the lowest two moments of the two-point function, 

S n,m,S:Cl,m',S'( E ^ W ) = p m,S^ p m',S'^ _ <Jp™,S(£)} |p-'^(^)| . (16) 

For (m, 5) = (to', S") they will give information about fluctuations and in particular about 
level motion in the ensemble [28]. For (m, S) ^ (to/, S'), the covariances (cross correlations) 
are non-zero for BEGOE while they will be zero for independent GOE representation for 
the to boson Hamiltonian matrices with different m or S. Note that the Q value has to be 
same for both (to, S) and (to/, S') systems so that the Hamiltonian in two-particle spaces 
remains same. Now we will discuss analytical and numerical results for S n and numerical 
results for S 22 for large values of (Q, m) and they are obtained using the results in Section 
IV. 

Trivially, the ensemble average of the energy centroids E c (m, S) will be zero [note that 
H is two-body for BEGOE(2)-s], i.e. (H) m ' s = 0. However the covariances in the energy 
centroids of H are non-zero and Eq. (8) gives, 



{H) m > b (H) r 

\2 \2 (17) 

Equations (13), (14) and (17) allow us to calculate S n for any (Q,m,S). For to = to' and 
S = S', the [S n] 1 ' 2 gives the width AE C of the fluctuations in the energy centroids. In the 
numerical calculations, we use Ao = Ai = A and therefore, Sn and S 22 are independent of A. 
Figure 6 gives some numerical results for AE C and it is seen that : (i) for to >> f2, the AE C 
is ~ 20% for S = and it goes down to ~ 15% for S = S max = to/2 for VI = 12; (ii) going 
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from ft = 12 to 40, AE C decreases to ~ 2 — 7%; (iii) for fixed (m, ft), there is decrease in AE C 
with increasing S value; (iv) for fixed (m, S) and very large m value, there is a sharp decrease 
in AE C with increasing ft up to ft ~ 20 and then it slowly converges to zero. It is possible 
to understand these results and the results for cross correlations [E n (m, S : m', S")] 1 / 2 , with 
(m, S) 7^ (to', S') as shown in Fig. 7, using the asymptotic structure of Q(Q,m, S). 

Let us consider the dense limit defined by m — > oo, ft — > oo and m/ft — > oo. Firstly the 
P u ' s {m, S) in Eq. (10) take the simpler forms, with S 2 = S{S + 1), 



™ =M=0 = K - 4£ 2 ) 2 p , =M=1 = 64m 2 ^ 2 (3m 2 + AS 2 ) 2 
32ft 4 ' 32ft 4 

=0 = (m 2 - 4^ 2 ) 2 =1 = 3m 4 + 40m 2 g 2 + 48(g 2 ) 2 

32ft 2 ' 32ft 2 

Using these in Eq. (13), with A = Ai = A, we have 

(m 2 + AS 2 ) 2 



(18) 



a 2 H (m, S) = A 2 



o 2 n [m } S)/ a 2 H (m, 
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m/(m + 2) + S 2 /S 2 



max 



(19) 



m/(m + 2) + 1 

The dense limit result given by Eq. (19) with m = 2000 is compared with the exact results 
in Fig. 5. Firstly, it should be noted that for the applicability of Eq. (19), ft should be 
sufficiently large and m » ft. Also, the result is independent of ft. Comparing with the 
ft = 12 and ft = 40 results, it is seen that the dense limit result is very close to the ft = 40 
results for m ^ 200. Thus for sufficiently large value of ft and m ^ 5ft, the dense limit 
result describes quite well the exact results. 



Simplifying (H) m,s (H) m ' s gives in the dilute limit, 



\2 20 
^ [{m 2 - AS 2 ) {(m>) 2 - A(S') 2 } + (3m 2 + AS 2 ) {3(m') 2 + A(S') 2 }] ■ 

Then pn] 1 / 2 , with m = rri and S = S' (for A = A x ) giving AE C , is 

ly u/2 a J? V2(5m* + 8m 2 S 2 + MS 2 ) 2 ) 

=AEc = ft (m 2 + AS 2 ) • (21) 

Eq. (21) gives [En] 1 / 2 to be y/lO/fl and 2/ft for S = and S = S max and these dense 
limit results are well verified by the results in Fig. 6b. Similarly, Eqs. (19) and (20) will 
give [En] 1/2 to be v^/ft for (m = rri : S = S max ,S' = 0) and 2/ft for (m = rri : S = 
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Smax, S' = S max — 1)- The upper and lower dashed lines in Fig. 7a for SI — 12 (similarly for 
S7 = 40) correspond to these two dense limit results respectively. It is seen that the dense 
limit results are close to exact results for SI = 40 but there are deviations for S7 = 12. Also, 
for fl = 40, the agreements are good only for m ^ 80 and these are similar to the results 
discussed earlier with reference to Fig. 5. 



Unlike for the covariances in energy centroids, we do not have at present complete ana- 
lytical formulation for the covariances in spectral variances. However, for a given mem- 
ber of BEGOE(2)-s, generating numerically (on a computer) the ensembles {V s=0 (2)} 
and {V s=1 (2)} and applying Eqs. (8) and (9) to each member of the ensemble will give 
{H 2 ) m ' s = a 2 (m,S) + [E c (m, S)} 2 . This procedure has been used with 500 members and 
results for E 22 are obtained for various (SI, m, S) values. For some examples, results are 
shown in Fig. 8 for both self correlations giving the width A (H 2 ) m " b of variances and cross 
correlations [E22] 1 / 2 with (m, S) 7^ (m',S'). It is seen that [E22] 1 / 2 are always much smaller 
than [En] 1 / 2 just as for EGOE(2) for spinless fermion systems [17]. It is seen from Fig. 8a 
that for il = 12, width of the fluctuations in the variances (H 2 ) m,S are ~ 3 — 5%. Simi- 
larly for large m, with Q very small, the widths are quite large but they decrease fast with 
increasing SI as seen from Fig. 8b. Finally, for S7 = 12, the cross correlations are ~ 4%. 



Besides the moments En and £22, it is possible to numerically construct the two-point 
function S n > m > S:a > m '> s ' {E, W) using the eigenvalues from the BEGOE(2)-s Hamiltonian ma- 
trices in small examples. We have carried out the calculations for S 4: ' 10 ' ' A ' 10 ' 1 (E,W) using 
a 500 member BEGOE(2)-s ensemble. It is seen that the structure of S(E, W) is similar to 
the nuclear shell model examples reported in [52] and EGOE(l+2)-s example in [22]. The 
maximum value of £4,10,0:4,10,1^ W ) is found to be _ 7% of p io,o(£) x p io,i(^)_ Let us add 

that it is important to identify measures involving En and E 22 and also S a ^ S:n ^ s ' (E , W), 
(m, S) 7^ (m', S') that can be tested using some experiments so that evidence for BEGOE(2) 
operation in real quantum systems can be established. 
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V. PREPONDERANCE OF S max = m/2 GROUND STATES AND NATURAL SPIN 
ORDER : ROLE OF EXCHANGE INTERACTION 

A. Introduction to regular structures with random interactions 

Johnson et al [53] discovered in 1998 that the nuclear shell model with random interactions 
generates, with high probability, + ground states in even-even nuclei (also generates odd- 
even staggering in binding energies, the seniority pairing gap etc.) and similarly, Bijker and 
Frank [40] found that the interacting boson model (sdlBM) of atomic nuclei [in this model, 
one considers identical bosons carrying angular momentum £ = (called s bosons) and £ = 2 
(called d bosons)] with random interactions generates vibrational and rotational structures 
with high probability. Starting with these, there are now many studies on regular structures 
in many-body systems generated by random interactions. See for example [54-56] for reviews 
on the subject. More recently, the effect of random interactions in the pn-sdlBM with F-spin 
quantum number has been studied by Yoshida et al [42]. Here, proton and neutron bosons 
are treated as the two components of a spin \ boson and this spin is called F-spin. Yoshida 
et al found that random interactions conserving F-spin generate predominance of maximum 
F-spin (F max ) ground states. It should be noted that the low-lying states generated by 
pn-sdTBM correspond to those of sdlBM and all scilBM states will have F = F max . Thus 
random interactions preserve the property that the low-lying states generated by pn-sdlBM 
are those of sdlBM. Similarly, using shell model with isospin conserving interactions (here 
protons and neutrons correspond to the two projections of isospin t = |), Kirson and 
Mizrahi [57] showed that random interactions generate natural isospin ordering. Denoting 
the lowest energy state (les) for a given many nucleon isospin T by Ei es (T), the natural 
isospin ordering corresponds to Ei es {T min ) < Ei es (T min + 1) < . . .; for even-even N=Z nuclei, 
T m in = 0. Therefore, one can ask if BEGOE(l+2)-s generates a spin ordering. 

As an application of BEGOE(l+2)-s, we present here results for the probability of gs 
spin to be S = S max and also for natural spin ordering (NSO). Here NSO corresponds to 
Eies(Smax) < Ei es (S max — 1) . . .. In this analysis, we add the Majorana force or the space 
exchange operator to the Hamiltonian in Eq. (1). Note that S in BEGOE(l+2)-s is similar 
to F-spin in the pn-sdlBM. First we will derive the exchange interaction and then present 
some numerical results. 
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B. U(fl) algebra and space exchange operator 



In terms of boson creation and annihilation (b) operators, the sp states for (f2) m 
systems are \i, m s ± |) = &t i m |0) with % — 1, 2, . . . , Q. It can be easily identified that the 
Ail 2 number of one-body operators A\j , 

^=( 6 ^X ; r=0,1 ' (22) 

generate C/(2Q) algebra. In Eq. (22), &.i ms = (-l)2 +T "-6 i) i _ ms . The U{2Q) irreducible 
representations (irreps) are denoted trivially by the particle number m as they must be 
symmetric irreps {m}. The Q 2 number of operators A ^ generate U(fl) algebra and similarly 
there is a U (2) algebra generated by the number operator h and the spin generators Sp, 

n = V2j2<; Sl = ^=J2 A ^- (23) 

i v i 

Then we have the group-subgroup algebra U (2f2) D U (Q) <S> SU (2) with SU(2) generated 
by As the U{2) irreps are two-rowed, the U(Q) irreps have to be two-rowed and they 
are labeled by {mi,m 2 } with m = mi + m 2 and S = (m 1 — m 2 )/2; mi > m 2 > 0. Thus 
with respect to U(Q) <g> SU(2) algebra, many boson states are labeled by |{m 1 ,m 2 },^) or 
equivalently by \{m, S),£), where £ are extra labels required for a complete specification of 
the states. The quadratic Casimir operator of the U (ft) algebra is, 

C 2 [U{n)]=2Y J A%-A% (24) 



ho 



and its eigenvalues are (C 2 [U(il)])^ mi ' m2 ^ = m 1 (m 1 +fi — l)+m 2 (m2 + fi- 3) or equivalently, 

(C 2 [U{Q)]) {m ' S) = ^(2fi + m - 4) + 25(5 + 1) . (25) 

Note that the Casimir invariant of SU{2) is S 2 with eigenvalues S(S+ 1). Now we will show 
that the space exchange or the Majorana operator M is simply related to C 2 \U(Q)\. 

Majorana operator M acting on a two-particle state exchanges the spatial coordinates of 
the particles (index i) and leaves the spin quantum numbers (m s ) unchanged. The operator 
form of M is 
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Equation (26) gives, with k a constant, 

M = | {C 2 [U(Q)\ - Qn} . (27) 

Then, combining Eqs. (25) and (27), we have 

M = K^n(j-lj + £ 2 j . (28) 

As seen from Eq. (28), exchange interaction with k > generates gs with S = S min = 0(|) 
for even(odd) m (this is opposite to the result for fermion systems where the exchange 
interaction generates gs with S = S max = m/2 [21, 22]). Now we will study the interplay 
between random interactions and the Majorana force in generating gs spin structure in boson 
systems. Note that for states with boson number fixed, M oc S 2 as seen from Eq. (28) and 
therefore, from now on, we refer to S 2 as the exchange interaction. 

C. Numerical results for S max = m/2 ground states and natural spin order 

In order to understand the gs structure in BEGOE(l+2)-s, we have studied P(S = S max ), 
the probability for the gs to be with spin S max = m/2, by adding the exchange term As S 2 
with \ s > to the Hamiltonian in Eq. (1), i.e. using 

{^} B EGOE(l+2)-S:Exch = Ml) + A [ {^=°(2)} + {V=\2)}] + \ s S 2 . (29) 

Note that the operator S 2 is simple in the (m, S) basis. Fig. 9a gives probability P(S = 
Smax) for the ground states to have spin S = S max as a function of exchange interaction 
strength \ s for A = Ai = A = 0, 0.1, 0.2, 0.3 and 0.5 and also for h(l) = with A = 1. 
Similarly, Fig. 9b shows the results for NSO. Calculations are carried out for (Q = 4, 
m = 10) system using a 500 member ensemble and the mean-field Hamiltonian h(l) is as 
defined in Section II. 

1. Preponderance of S max = m/2 ground states 

Let us begin with pure random two-body interactions. Then h(l) = in Eq. (29). 
Now in the absence of the exchange interaction (X s = 0), as seen from Fig. 9a, ground 
states will have S = S max , i.e. the probability P(S = S max ) = 1. The variance propagator 
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(see Fig. 5) derived earlier gives a simple explanation for this by applying the Jacquod 
and Stone prescription as discussed in Sec. IV B. Thus pure random interactions generate 
preponderance of S — S max ground states. On the other hand, as discussed in Section V B, 
the exchange interaction acts in opposite direction by generating S = S m i n ground states. 
Therefore, by adding the exchange interaction to the {l / (2)} ensemble, P(S = S, max ) starts 
decreasing as the strength As (As > 0) starts increasing. For the example considered in Fig. 
9a, for \ s > 4, we have P(S = S max ) ~ 0. The complete variation with \ s is shown in Fig. 
9a marked h(l) = and A = 1. 

Similarly, on the other end, for A = in Eq. (29), we have H — h(l) in the absence of 
the exchange interaction. In this situation, as all the bosons can occupy the lowest sp state, 
gs spin S = S m ax- Therefore, P(S = S max ) = 1. When the exchange interaction is turned 
on, P(S = S max ) remains unity until As equals the spacing between the lowest two sp states 
divided by to. As in our example, the sp energies are q = i + l/i, we have P(S = S max ) = 1 
for As < 0.05. Then P(S = S max ) drops to zero for As > 0.05. This variation with As is 
shown in Fig. 9a marked A = 0. Figure 9a also shows the variation of P(S = S max ) with 
As for several values of A between 0.1 and 0.5. It is seen that there is a critical value (X c s ) 
of As after which P(S = S max ) = and its value increases with A. Also, the variation of 
P(S = S max ) with As becomes slower as A increases. 

In summary, results in Fig. 9a clearly show that with random interactions there is 
preponderance of S = S max — to/2 ground states. This is unlike for fermions where there 
is preponderance of S — S m i n = 0(|) ground states for to even(odd). With the addition of 
the exchange interaction, P(S = S max ) decreases and finally goes to zero for As > A 5 and 
the value of \ c s increases with A. We have also carried out calculations for (Q — 4, to = 11) 
system using a 100 member ensemble and the results are close to those given in Fig. 9a. 
All these explain the results given in [42] where random interactions are employed within 
pn-sdlBM. 

2. Natural spin ordering 

For the system considered in Fig. 9a, for each member of the ensemble, eigenvalue of the 
lowest state for each spin S is calculated and using these, we have obtained total number of 
members N x having NSO as a function of As for A = 0.1, 0.2 and 0.3 using the Hamiltonian 
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given in Eq. (29). As stated in Section V A, the NSO here corresponds to (as S = S max is 
the spin of the gs of the system) Ei es (S max ) < Ei es (S max - 1) < E les (S max -2) < .... The 
probability for NSO is N x /500 and the results are shown in Fig. 9b. In the absence of the 
exchange interaction, as seen from Fig. 9b, NSO is found in all the members independent 
of A. Thus random interactions strongly favor NSO. The presence of exchange interaction 
reduces the probability for NSO. Comparing Figs. 9a and 9b, it is clearly seen that with 
increasing exchange interaction strength, probability for gs state spin to be S = S max is 
preserved for much larger values of As (with a fixed A) compared to the NSO. Therefore for 
preserving both S = S max gs and the NSO with high probability, the \ s value has to be 
small. We have also verified this for the (Q — 4, m — 11) system. Finally, it is plausible 
to argue that the results in Fig. 9 obtained using BEGOE(l+2)-s are generic for boson 
systems with spin. Now we will turn to pairing in BEGOE(2)-s. 

VI. PAIRING IN BEGOE(2)-s 

Pairing correlations are known to be important not only for fermion systems but also 
for boson systems [46]. An important issue that is raised in the recent years is: to what 
extent random interactions carry features of pairing. See [6, 54, 55, 58] for some results for 
fermion systems. In order to address this question for boson systems, first we will identify 
the pairing algebra in (fi, m, S) spaces of BEGOE(2)-s. Then we will consider expectation 
values of the pairing Hamiltonian in the eigenstates generated by BEGOE(2)-s as they carry 
signatures of pairing. 

A. U(2n) D [U(Q) D SO (SI)] SU S {2) Pairing symmetry 

In constructing BEGOE(2)-s, it is assumed that spin is a good symmetry and thus the 
m-particle states carry spin (S) quantum number. Now, following the 5*0(5) pairing algebra 
for fermions [59], it is possible to consider pairs that are vectors in spin space. The pair 
creation operators P i:M for the level i and the generalized pair creation operators (over the 
Q levels) P^, with /i = —1, 0, 1, in spin coupled representation, are 

p, = (MX = E^, • to* = ^E-D 1 -" (<*)'„ ■ po) 
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Therefore in the space defining BEGOE(2)-s, the pairing Hamiltonian H p and its two-particle 
matrix elements are, 

h p = Y. p » ( p rf ' ^ s i h p i fa» = ^ ^ 5 w • ( 31 ) 

With this, we will proceed to identify and analyze the pairing algebra. It is easy to verify 
that the Q(Q — 1) /2 number of operators CV,- = A^ — A ^, i > j generate a SO(Q) subalgebra 
of the U(Q) algebra. Therefore we have U(2Q) D [U(Q) D SO(tt)} <g> SU(2). We will show 
that the irreps of SO (CI) algebra are uniquely labeled by the seniority quantum number v 
and a reduced spin s similar to the reduced isospin introduced in the context of nuclear shell 
model [60] and they in turn define the eigenvalues of H p . The quadratic Casimir operator 
of the SO (CI) algebra is, 

C 2 [SO(a)] = 2j2Cv-C Ji . (32) 

i>j 

Carrying out angular momentum algebra [61] it can be shown that, 

C 2 [SO(Q)] = C 2 [U(Q)]-2H p -n . (33) 

The quadratic Casimir operator of the U(Q) algebra is given in Eq. (24). Before discussing 
the eigenvalues of the pairing Hamiltonian H pi let us first consider the irreps of SO(Ct). 

Given the two-rowed U(Q) irreps {mi,m 2 }; mi + m 2 = m, mi — m 2 = 2S, it should 
be clear that the SO(fl) irreps should be of [t>i,t>2] type and for later simplicity we use 
Vi + v 2 = v and v\ — v 2 — 2s. The quantum number v is called seniority and s is called 
reduced spin. The SO(Q) irreps for a given {m 1 ,m 2 } can be obtained as follows. First 
expand the U(Q) irrep {m 1 ,m 2 } in terms of totally symmetric irreps, 

{mi,m 2 } = {mi} x {m 2 } - {m 1 + 1} x {m 2 - 1} . (34) 

Note that the irrep multiplication in Eq. (34) is a Kronecker multiplication [62, 63]. For a 
totally symmetric U(Q) irrep {m 1 }, the SO(Q) irreps are given by the well-known result 

{m} ->■ [v] = [m] © [rri - 2] © ... © [0] or [1] . (35) 

Finally, reduction of the Kronecker product of two symmetric SO(Q) irreps [vi] and [v 2 ], 
Q > 3 into SO(Q) irreps [t>i,t>2] is given by (for vi > v 2 ) [62, 63], 

1)2 v 2 —k 

[vi] x [v 2 \ = ^^[vi-v 2 + k + 2r } k}@ . (36) 

k=0 r=0 
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Combining Eqs. (34), (35) and (36) gives the {mi,m 2 } — > [vi,v 2 ] reductions. It is easy to 
implement this procedure on a computer. 

Given the space defined by \{mi, m2}, [i>i, 1*2], ex), with a denoting extra labels needed for 
a complete specification of the state, the eigenvalues of C 2 [SO( c T)] are [62] 

(C 2 [SO(Q)]) {mi ' m2} ' [vi ' V2] =v 1 (v 1 + Q-2) + v 2 (v 2 + Q - 4) . (37) 

Now changing {mi,m 2 } to (m,S) and [t>i,t>2] to (v,s) and using Eqs. (33) and (25) will 
give the formula for the eigenvalues of the pairing Hamiltonian H p . The final result is, 

E p (m,S,v,s) = (H p ) m ' s ' v '~ s = j(m-v)(2n-6 + m + v) + [S(S + l)-s(s + l)] . (38) 

This is same as the result that follows from Eq. (18) of [59] for fermions by using Q — > —Q 
symmetry. From now on, we denote the U(Q) irreps by (to, 5) and SO (ft) irreps by (v, s). 
In Table I, for (Q,m) = (4, 10), (5, 8) and (6,6) systems, given are the (m,S) — > (v,s) 
reductions, the pairing eigenvalues given by Eq. (38) in the spaces defined by these irreps 
and also the dimensions of the U(Q) and SO(Q) irreps. The dimensions d(Q,,m, S) of the 
U(Vt) irreps (m,S) are given by Eq. (2). Similarly, the dimension d(vi,v 2 ) d(v,s) of the 
SO(Q) irreps [^1,^2] follow from Eqs. (35) and (36) and they will give 

V2 — 1 V2~k 

d(vi,v 2 ) = d(vi)d(v 2 ) - ^ ^ d(vi - v 2 + k + 2r, k) ; 

fc=0 r=0 

(39) 



d(v) = 



v J V v - 2 



Note that in general the SO (ft) irreps (v, s) can appear more than once in the reduction of 
U(Cl) irreps (m, S). For example, (2, 1) irrep of SO(tt) appears twice in the reduction of the 
U{Q) irrep (10,1). 

It is useful to remark that just as the fermionic 5*0(5) pairing algebra for nucleons in j 
orbits [59, 64, 65], there will be a 5*0(4, 1) complementary pairing algebra corresponding to 
the SO(Q) subalgebra. The ten operators P^, (P^Y, S^ and n form the 50(4, 1) algebra. 
It is possible to exploit this algebra to derive properties of the eigenstates defined by the 
pairing Hamiltonian but this will be discussed elsewhere. 
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B. Pairing expectation values 



Pairing expectation values are defined by (H P ) S ' E = (to, S, E \ H p | to, S, E) for eigen- 
states with energy E and spin S generated by a Hamiltonian H for a system of to bosons in 

S E 

fl number of sp orbitals (for simplicity, we have dropped Q and to labels in (H p ) ' ). In our 
analysis, H is a member of BEG0E(2)-s. As we will be comparing the results for all spins 
at a given energy E, for each member of the ensemble the eigenvalues for all spins are zero 
centered and normalized using the m-particle energy centroid E c {m) = (H) m and spectrum 
width a(m) = [{H 2 ) m — {^(m)} 2 ] 1 / 2 . Then the eigenvalues E for all S are changed to 
E — [E — E c (m)]/a(m). Using the method described in Section II, the H p matrix is con- 
structed in good Ms basis and transformed into the eigenbasis of a given S for each member 
of the BEG0E(2)-s ensemble. Then the ensemble average of the diagonal elements of the 

S E S £/ 

H p matrix will give the ensemble averaged pairing expectation values (H p ) ' <^ (H p ) ' . 

Using this procedure for a 500 member BEGOE(2)-s ensemble with = 4, to — 10 and 

s E ^ ^ 

S = — 5, results for (H p ) ' as a function of energy E (with E as described above) and 

spin S are obtained and they are shown as a 3D histogram in Fig. 10. From Table I, it is 

seen that the maximum value of the eigenvalues E p (m, S,v, s) increases with spin S for a 

&xed-(n,m). The values are 28, 32, 34, 42, 48 and 60 for S = - 5 respectively for fl = 4 

and m = 10. Numerical results in Fig. 10 also show that for states near the lowest E 
s E 

value, (H p ) ' increases with spin S. Thus random interactions preserve this property of 
the pairing Hamiltonian in addition to generating S = S max ground states as discussed in 
Section VC. It is useful to remark that random interactions will not generate S = S max 
ground states with (v, s) = (m, m/2) as required for example in the pn-sdlBM. This needs 
explicit inclusion of pairing and exchange terms in the Hamiltonians defined by Eqs. (1) 
and (6). 

For a given spin S, the pairing expectation values as a function of E are expected, for 
two-body ensembles, to be given by a ratio of expectation value density (EVD) Gaussian (the 
first two moments given by (H p H} m,s and (H p H 2 } m ' S ) and the eigenvalue density Gaussian 
with normalization given by (H p ) m ' s and this itself will be a Gaussian [18]. Let us denote 
the EVD centroid by E c (m, S : H p ) and width by a(m, S : H p ). Then the ratio of Gaussians 
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will give 




S,E 




cxp 



2[1 -a 2 (m, S)] 



e*(m, S) 



cxp 



Q\m,S)-l) 
2d 2 {m,S) 



E — 



1 — a 2 (m, S) 



?(m, S) 




(40) 



Here, e(m, S) = {E c (m, S : H p ) — E c {m, S)}/a(m, S), a(m, S) = a(m, S : H p )/a(m, S) and 
E = [a(m)/a(m, S)]{E — £}; 8 = [E c (m,S) — E c (m)]/a(m). The Gaussian form given 
by Eq. (40) is clearly seen in Fig. 10 and this also gives a quantitative description of the 
results. Note that in our example, e{10,S) = 0.001, 0.001, 0.001, 0.002, 0.002, 0.003 and 
£(10, S) = 1.045, 1.047, 1.053, 1.062, 1.073, 1.082 respectively for S = - 5. 

VII. CONCLUSIONS 

In the present work, we have introduced the BEGOE(l+2)-s ensemble and a method 
for constructing BEGOE(l+2)-s for numerical calculations has been described. Numerical 
examples are used to show that, like the spinless BEGOE(l+2), the spin BEGOE(l+2)-s 
ensemble also generates Gaussian density of states in the dense limit. Similarly, BEGOE(2)- 
s exhibits GOE level fluctuations. On the other hand, BEGOE(l+2)-s exhibits Poisson 
to GOE transition as the interaction strength A is increased and the transition marker 
A c is found to decrease with increasing spin. Moreover, ensemble averaged covariances 
in energy centroids and spectral variances for BEGOE(2)-s between spectra with different 
particle numbers and spins are studied using the propagation formulas derived for the energy 
centroids and spectral variances. For Q = 12 systems, the cross correlations in energy 
centroids are ~ 15% and they reduce to ~ 4% for spectral variances. We have also derived the 
exact formula for the ensemble averaged fixed- (m, S) spectral variances and demonstrated 
that the variance propagator gives a simple explanation for the preponderance of spin S = 
Smax ground states generated by random interactions as in pn-sdlBM. It is also shown, by 
including exchange interaction S 2 in BEGOE(l+2)-s, that random interactions preserving 
spin symmetry strongly favor NSO (just as with isospin in nuclear shell model). In addition, 
we have identified the pairing SO(Q) symmetry and showed using numerical examples that 
random interactions exhibit pairing correlations in the gs region and also they generate a 
Gaussian form for the variation of the pairing expectation values with respect to energy. 

Results in this paper represent a beginning in systematic studies of BEGOE(l+2)-s 
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ensemble. To go beyond the present work requires: (i) dealing with systems with much 
larger Q and to values and then the matrix dimensions will be 10 5 — 10 6 and higher; (ii) 
further group theoretical analysis for deriving higher moments of the two-point function and 
this requires new results for Wigner and Racah coefficients for general U (Q) algebras (see for 
example [11]); and (iii) identifying and analyzing measures that can be used in experimental 
data analysis, say for example data from ultracold gases. With progress in (i) and (ii), it 
is possible to investigate more systematically BEGOE(l+2)-s as a function of the strength 
parameter A by analyzing the derealization and other related measures (see for example 
[6, 25, 26]). With these it is possible to address (iii) but this is for future. 

Further extensions of BEGOE(l+2)-s including s = 1, 2, . . . degrees of freedom for 
bosons, as emphasized in the Introduction, are relevant for spinor BEC studies [46, 47]. 
These extended BEGOE's will be explored in future. Finally, it is useful to mention that 
the numerical code developed for constructing BEGOE(l+2)-s ensemble in fixed- (to, S) 
spaces can be used in analyzing ensembles generated by Hamiltonians preserving only Ms- 
Here it is important to point out that, besides the SO(Q) pairing discussed in Section VI, 
there is another pairing algebra that corresponds to U(2Q) D SO(2Q) D [SO(Q) <g> 5*0(2)] 
as pointed out in [62] . The pairing Hamiltonian corresponding to this group-subgroup chain 
preserves only Ms but not S and hence, relevant for BEG0E(l+2) with fixed Ms- 
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APPENDIX A 

Let us consider a system of to spinless bosons occupying N sp states | i/j) , i — 1, 2, . . . , N 
and the Hamiltonian is say two-body. Then the Hamiltonian operator is 




(Al) 
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with the symmetries for the symmetrized two-body matrix elements (v k v\ \ H \ Vj) being, 

(v k vi\H\ Vj i>i) = {v k v x | H | vi v>) , 

{A2) 

(v k vi \ H | v,i Vj) = (vi Vj | H | v k vi) . 

The Hamiltonian matrix H(m) in m-particle spaces contains three different types of non-zero 
matrix elements and explicit formulas for these are 



( n (*r i h i n Kr) = E ni n' +s s f (vw i h i v t v 3 ) , 

\r=i,j,... r=i,j,... / i>j V %3 ' 

Uv t r-\vjf^ n {v,r\H\ n (v r r\ = 

\ r'=k,l,... r=i,j,... / 



rij (rij + 1) (% - 4'^)' 
(1 + <^) (1 + 8 V j) 



1/2 



(A3) 



\ r'=m,n,... r=i,j,... / 



n k (ni - 5 M ) (rii + 1) (n^ + 1 + (%) 



nl/2 



(l/jl/j | H | Z/fcZ/;) . 



(i + <y (1 + 4z) 

In the second equation in Eq. (A3), i ^ j and in the third equation, four combinations are 
possible: (i) k — I, i — j, k ^ i; (ii) k = I, i ^ j, k ^ i, k ^ j; (iii) k ^ I, i = j, i ^ k, 
i ^ I; and (iv) i ^ j ^ k ^ I. BEGOE(2) for spinless boson systems is defined by Eqs. (A2) 
and (A3) with the H matrix in two-particle spaces being GOE. Note that the H{m) matrix 
dimension is ^ N+m ~ 1> j and the number of independent matrix elements is ^2(^2 + 1)/2 where 
(I2 = N(N + l)/2. It is useful to mention that the formulas for the energy centroids and 
spectral variances in m-boson spaces follow from Eqs. (8), (9), (10) and (11) with S = m/2, 
putting the s = matrix elements to zero and replacing Q by N. 
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TABLE I: Classification of states in the U(2Q) D [U(tt) D SO(Q)]® SU S {2) limit for (0, m) = (4, 10), (5, 8) 
and (6,6). Given are U(Cl) labels (m,S) and SO(Q) labels (v,s) with the corresponding dimensions 
d(p., m, S) and d(v, s) respectively and also the pairing eigenvalues E p (m, S, v, s). Note that s rd(v, s) = 
d(Q,m,S); here r denotes multiplicity of the SO (ft) irreps and in the table, they are shown only for the 
cases when r > 1. 
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FIG. 1. (Color online) Ensemble averaged eigenvalue density p m,s (E) vs E for BEGOE(2)-s 
ensembles with = 4, m = 10 and O = 4, m = 11. In the figure, histograms constructed with 
a bin size 0.2 are BEGOE(2)-s results and they are compared with Gaussian (dashed red) and 
Edgeworth (ED) corrected Gaussian (solid green) forms. The ensemble averaged values of the 
excess parameter (72) are also shown in the figure. In the plots, the area under the curves is 
normalized to the dimensions d(Q,m, S). See text for further details. 
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FIG. 2. (Color online) (a) Ensemble averaged nearest neighbor spacing distribution (NNSD) and 
(b) Dyson- Mehta statistic A^(L) vs L for L < 60. Results are for the same systems considered 
in Fig. 1; first column gives the results for (f2 = 4, m = 10) and the second column for (O = 4, 
m = 11) systems. The NNSD histograms from BEGOE(2)-s are compared with Poisson (dashed 
red) and GOE (Wigner) forms (solid green) and similarly the A^(L) results. In the NNSD graphs, 
the bin-size is 0.2 and x is the nearest neighbor spacing in the units of local mean spacing. See 
text and Fig. 1 for further details. 



36 



to 
c 

<D 

T5 

<D 
3 

> 
c 

<D 

5> 
Li 



(a) 



>< 

ST 



(b) 



BEGOE(1+2)-s; X=0A 
Q.=4, m=1 1 , S=7/2 




0.8 
0.6 
0.4 
0.2 
0.0 



Poisson 
GOE 



(c) 



20 40 
L 



60 



FIG. 3. (Color online) (a) Ensemble averag ed eigenvalue density p m ' S {E), (b) NNSD and (c) A 3 (L) 
vs. L for a 100 member BEGOE(l+2)-s ensemble for Q = 4, m = 11 and S 1 = 7/2 system with 
Aq = Ai = A = 0.1 in Eq. (1). For all other details, see text and Figs. 1 and 2. 
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FIG. 4. (Color online) NNSD for a 100 member BEGOE(l+2)-s ensemble with Q = 4, m = 10 
and spins S = 0, 2 and 5. Calculated NNSD are compared to the Poisson (red dashed) and Wigner 
(GOE) (green solid) forms. Values of the interaction strength A and the transition parameter 
A are given in the figure. The values of A are deduced as discussed in [6]. The chaos marker A c 
corresponds to A = 0.3 and its values, as shown in the figure, are 0.039, 0.0315, 0.0275 for S = 0, 2 
and 5 respectively. Bin-size for the histograms is 0.2. 
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FIG. 5. (Color online) BEGOE(2)-s variance propagator Q(fi,m,5)/Q(fi, rn,S max ) vs S/S max for 
various values of f2 and m. Formula for Q(f2, m,S) follows from Eqs. (10), (13) and (14). Note 
that the results in the figure are for Ao = Ai = A in Eq. (6) and therefore independent of A. Dense 
limit (dot-dashed) curve corresponds to the result due to Eq. (19) with m = 2000. 
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FIG. 6. (Color online) (a) Self-correlations E-q in energy centroids, giving width AE C of the 
fluctuations in energy centroids scaled to the spectrum width, as a function of spin S for different 
values of m and f2. (b) Self correlations as a function of f2 for 5000 bosons with minimum spin 
(S = 0) and maximum spin (S = 2500). Dense limit (dot-dashed) curves for 5 = and S = m/2 
in (b) correspond to the results given by Eq. (21). See text for details. 
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FIG. 7. (Color online) Cross-correlations in energy centroids for various BEGOE(2)-s systems, 
(a) S^{ 2 vs m with m = mf but different spins (5 7^ S"). (b) S^{ 2 vs m! with m = 10 and S = S' = 
and S = 5, S' = m'/2. The dashed lines in (a) are the dense limit results. See text for details. 
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FIG. 8. (Color online) Correlations in spectral variances £ 2 2 I * or var i° us BEGOE(2)-s systems, 
(a) Self-correlations, giving width A (H 2 ^ m,S of the spectral variances, as a function of spin S 
for to = 12, 24 and 120 with = 12. (b) Self-correlations as a function of Q for 5000 bosons 
with S = and 2500. (c) Three examples for cross-correlation in spectral variances with same 
or different particle numbers and same or different spins. All the results are obtained using 500 
member ensembles. See text for details. 
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FIG. 9. (Color online) (a) Probability for ground states to have spin 5 = S m ax as a function of the 
exchange interaction strength As > 0. (b) Probability for natural spin order (NSO) as a function 
of As- Results are shown for a 500 member BEGOE(l+2)-s : Exch ensemble generated by Eq. 
(29) for a system with = 4 and m = 10. Values of the interaction strength A are shown in the 
figure. 
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FIG. 10. (Color online) Ensemble averaged pairing expectation values (H p ) ' vs E and S, shown 
as a 3D histogram, for a 500 member BEGOE(2)-s ensemble with f2 = 4 and m = 10. The bin-size 
is 0.2 for E. Note that the E label in this figure is different from the E used in Figs. 1 and 3a. 
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